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Abstract 

The behaviour of large amphtude electrostatic waves in cold plasma is investigated in the context of Born- 
Infeld electrodynamics. The equations of motion for the plasma are established using an unconstrained action 
principle. The maximum amplitude and frequency of a large amplitude electrostatic wave are determined, 
and a lower bound on the wavelength is established. The maximum electric field is found to be the same as 
that on a point electron at rest. 

1 Introduction 

Recent years have seen renewed interest in non-linear electrodynamics [IHS] , in particular Born-Infeld theory |3] , 
and its implications. Unlike their Maxwell counterparts, the Born-Infeld field equations are fundamentally non- 
linear; they are non-linear even in the vacuum. Born and Infeld introduced their theory [J in order to ameliorate 
the singular self-energy of a point charge. Furthermore, it was discovered that among the family of non-linear 
generalizations of Maxwell electrodynamics, Born-Infeld theory possesses a number of highly attractive features; 
in particular, like the vacuum Maxwell equations, the vacuum Born-Infeld equations exhibit zero birefringence 
and its solutions have exceptional causal behaviour [SJIH]. Furthermore, Born-Infeld theory shares a number of 
properties with the low energy dynamics of strings and branes [7] . 

Some of the most extreme conditions ever encountered in a terrestrial laboratory are created when high- 
power laser pulses interact with matter. The laser pulse immediately vaporizes the matter to form an intense 
laser-plasma providing novel avenues for generating intense bursts of coherent electromagnetic radiation for a 
wide range of applications in biological and material science [5]. Furthermore, laser-plasmas permit controllable 
investigation of matter in extreme conditions that only occur naturally away from the Earth. It is expected 
that the next generation of ultra-intense lasers will, for the first time, allow controllable access to regimes 
where a host of different quantum clcctrodynamic phenomena will be evident and theories of non-linear 
electrodynamics will be central to their study. However, it is conceivable that Born-Infeld electrodynamics may 
play a role before it is necessary to invoke quantum electrodynamics [3]. 

Motivated by recent analyses of the propagation of Born-Infeld electromagnetic waves in waveguides [5] and 
in background uniform magnetic fields |3] , this article is addresses aspects of the non-linear electrodynamics of 
laser-plasmas. 

A sufficiently short and intense laser pulse propagating through a plasma may create a travelling longitudinal 
plasma wave whose velocity is approximately the same as the laser pulse's group velocity. However, it is not 
possible to sustain arbitrarily large electric fields; substantial numbers of plasma electrons become trapped in 
the wave and are accelerated, which dampens the wave (the wave 'breaks'). Early theoretical investigation of 
non-linear plasma waves was undertaken by Akhiezer and Polovin [10| , and later expounded by Dawson in 
the context of wave-breaking. 

Wave-breaking is a fundamentally non-linear phenomenon, and it is natural to explore the properties of 
Born-Infeld plasmas from this perspective. For simplicity, in the present article we adopt the cold plasma model 
for the electron fluid and generalize well-known results [TUIIII], for the cold Maxwell plasma, to the context 
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of non-linear electrodynamics. In particular, we obtain an exact expression for the maximum electric field 
of electrostatic oscillations in a cold Born-Infeld plasma and analyse the wavelength and frequency of those 
oscillations. There are numerous discussions of electrostatic waves in a warm Maxwell plasma; for example, 
see [I1IT3]. 

We employ the Einstein summation convention throughout this article. Latin indices a, &, c run over 0, 1, 2, 3 
and units are used in which the speed of light c = 1 and the permittivity of the vacuum Eq ^ 1. 

Let (x") be an incrtial coordinate system on Minkowski spacctimc {Ai,g) where a;*^ is the proper time of 
observers at fixed Cartesian coordinates {x^,x'^,x^) in the laboratory. The metric tensor g has the form 

g = i^abdx" ®dx^ (1) 

with 

-lifa = & = 0, 

Tlab^ { 1 if a = 5 ^ 0, (2) 

and the Hodge map ★ is induced from the 4-form ^1 on Ai where 

-kl = dx^ A dx^ ^ dx^ ^ dx^ . (3) 

The plasma electrons are represented as a cold relativistic fluid; their worldlines are trajectories of the unit 
normalized future-pointing timelike 4-vector field V oti M and the 0-form n is their proper number density. 
We are interested in the evolution of a plasma over timescales during which the motion of the ions is negligible 
in comparison with the motion of the electrons, and we assume that the ions are at rest and distributed 
homogeneously in the laboratory frame. Their worldlines are trajectories of the vector field iVjon = niond/dx^ 
on M. where riion is the ion number density (a positive constant) in the laboratory frame. 

Maxwell's equations may be written covariantly as 

dF = 0, d-kG ^ -qn *V - qion * Mon (4) 

where g < is the charge on the electron and gion = Z\q\ is the charge on an ion with Z the multiplicity of 
the ionization. The 1-forms V and iVjon are the metric duals of the vector fields V and A^ion respectively, i.e. 
the 1-form V satisfies V{U) = g{V, U) for all vector fields U on M. The Maxwell 2-form F encodes the electric 
field E and the magnetic induction B; the excitation 2-form G encodes the electric displacement D and the 
magnetic field H. In Maxwell electrodynamics G = -F in the vacuum; however, in Born-Infeld electrodynamics 
(m is retained but the constitutive relations for (D, H) in terms of (E, B) arc non-linear in the vacuum. In 
Born-Infeld electrodynamics the constitutive relation for G in terms of F is 

G= , ^ (F-k'^Y*F) (5) 

where the invariants X and Y are 

X = k{FA*F), Y = k{FAF) (6) 

and At is a new constant of nature. 

To motivate the field equations used herein, we begin section [2] by developing them using an unconstrained 
action principle [T3]. We will briefly review how ([5]) arises and show that V satisfies the field equations 

^yV=±,yF, g{V,V)^-l (7) 
m 

for a relativistic cold electron fluid, where qtyF is the Lorentz 4-force acting on the electron fluid, by is the 
interior product with respect to V , m is the electron rest mass and V is the Levi-Civita connection on Ai. 
This will be followed in section |2] by an exploration of the properties of large amplitude relativistic electrostatic 
oscillations. 
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2 Action principle 



A succinct action principle for establishing (jH [7]) is formulated using a i/iree-dimensional manifold B (often 
called a 'material' jl5j . or 'body', manifold), where each point in B corresponds to an integral curve of V. 

The independent variables of the action introduced below are the electromagnetic 1-form A on A4 and a 
submersion / from Ai to B, 

f:M ^ B (8) 
x'^ ^ ^^^f\x). (9) 

The map / identifies an integral curve of V with its corresponding point in B, where upper case Latin indices 
j4, i? = 1, 2, 3. The inverse image f~^(P) of the point P G S is the worldlinc of an idealized 'material' particle in 
the electron fluid and V may be determined from / as the unique timelikc future-pointing vector field satisfying 

dJ^{V) = (10) 

and 

g{V,V)^-l. (11) 

In the present approach, material properties of the electron fluid are encoded in tensors and forms on S, 
and fields on M. are induced from material tensors on B using /. In particular, the number of electrons -/V[E] 
occupying the spacelike hypersurfacc S may be induced from a 3-form f7 on B as 

ivp] = / (12) 

where /* denotes the pull-back map, induced from /, on forms on B to forms on lA. The electron number 
current 3-form j , 

3 = (13) 

is identically closed, 

dj = 0, (14) 

where d is the exterior derivative, which follows because is a top form on B and df* = f*d. 

A variational principle yielding ([7]) for V and (jH) for the potential 1-form A, where F = dA, employs the 
action functional S where 

S[AJ]= [ [CEM{X,Y)*1^7n^r]*l~qA/\j^qionAAjion)]- (15) 

The 0-form CEuiX, Y) is a local function of the invariants X and Y only. For standard Maxwell electrodynamics 
^EM ~ -^/2, whereas for Born-Infcld electrodynamics £em = 'Cbi and 
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CbiIX, Y) = —(1 - ^l-K^X^K'^Y^A). (16) 



The 3-form jion = nmndx^ A dx^ A dx^ is the background ion number current and the 0-form j ■ j is the square 
of the magnitude of j , 

j -J =*~\j A^j) (17) 

where is the inverse of the Hodge map 

The first term £em(-^, F) ★ 1 in the Lagrangian 4-form in (ITSi) depends only on A and g, the second term 
m^/j ■ j -k 1 depends only on / and g, the third term qA A j couples / and A and the fourth term qionA A jion 
depends only on A. 
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We will show in the following that the field equations (jH [7]) for n, V and A are recovered from those for / 
and A by introducing the electron proper number density n and 4- velocity field V as 



V 



(18) 
(19) 



2.1 Non-linear generalization of the Maxwell equations 

Equations for the electromagnetic field F arise upon seeking stationary variations of S with respect to A 



SaS^ 



M 



■^aX + SaY ] -kl - 5AA {qj + QionJion) 



V dX 



dY 



where Saoc denotes the variation of a form a with respect to A. Using 1^ and -k-ka 
on A^, it follows 



X kl = -F AkF, 



Y*l = -FAF 



and 



SaX * 1 = 5a{X ★ 1) = -2d6A A kF, 6aY ★ 1 = 6a{Y * 1) = -2dSA A F. 



Hence 



5aS^ 



[ 6AA 




Im 





■*F -I ^T^F ] - [qj + (7ionJion) 



dX 



dY 



(20) 

-a where a is a 4-form 
(21) 

(22) 
(23) 



where the variation SA is chosen to have compact support on A4 , and an integration by parts has been used in 
the final step. By demanding that the action is stationary under all such variations we obtain 



dkG = -qj - qionju 



where 



-k F -\ F 



dX 



dY 



(24) 



(25) 



The field equations @ are obtained by introducing n and V using ([TOI and , and the ion number 4-current 

TVion as 



(26) 



and noting dF = since F is an exact 2-form. Equation ^ immediately follows from (fTHl [25)) with £em — )^b\ 
and using k:k: /3 ~ — /? where /3 is a 2-form. 



2.2 Field equations for the electron fluid 

The field equations ^ describing the electron fluid arc obtained by seeking stationary variations of S with 
respect to /. Using (fT8|) it follows 



'] - J 



-'{SfjAkj)] 



To proceed further we express j = f*n explicitly in terms of the components {/^} of /. 



1 /an 



ABC 



3!V 

Lwdj + diwj 



o f ] Sf"" df^ A df /\df' + - {^ABC o f) dSf" A df"" A df 



(27) 



(28) 
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where n = l^ABciO A d^^ A d^'^ and 



j = f*^ = ^ {nABC o /) df^ A rf/^ A d/^ 
have been used. The vector field W is 
W = &!^Wa 



(29) 



(30) 



where {H^a} is a basis for ^-orthogonal vector fields on A^, and the frame {V, Wa\ and coframe {— d/"*} are 
naturally dual, i.e. 



d!^{WB)=5i, Y{Wa)^^, df^iV)^0, ViV) = -l 
where Sg is the Kronecker delta. Since j is a closed 3-forni it follows 

5fj = dbwj 
using (pSj) and so 



diwj A 



1 



(31) 
(32) 

(33) 



Varying ([T5]) with respect to / yields 
SfS^ 



M 



j=jSfj A ★j + qAA Sfj 



where the variation 6f^ is chosen to have compact support on A4. An integration by parts yields 



SfS^ Lwj Ad{^=*j - qA) 
Jm V j ■ .7 

and by demanding that S is stationary under all suitable variations 5f^ we obtain 



i-WaJ a d{- 



/] -3 



*j-qA)=0. 



Eliminating j in favour of n and V using (|18l I19|) yields 
LwA^'vi'iTT'dV — qF) = 
and since {V, Wa} is a frame on M it follows 

LydV = —LyF. 



(34) 



(35) 



(36) 



(37) 



(38) 



Since g{V,V) is constant it may be shown LydV = S/yV (see, for example, [11]) and ([7]) follows immediately. 

3 Non-linear electrostatic oscillations 

The remainder of this article focuses on properties of large-amplitude longitudinal electrostatic waves propa- 
gating parallel to the x'^-axis with phase velocity v (with < w < 1) in the laboratory frame. We introduce the 
pair {e^, e^} 



;^ = vdx^ - dx° 



dx^ - vdx° 



(39) 



and note that the orthonormal coframe {76^, 76^, dx^, dx^} is adapted to observers moving at velocity v along 
(i.e observers in the 'wave frame'). We seek a 4- velocity V of the form 



y = M(C)e' + V'(C)f 



(40) 



where C = x"^ — fx" is the wave's phase and we have adopted the so-called 'quasi-static approximation'; the 
pointwise dependence of /i and ijj is on ( only. Using g{V, V) = —I the component ip is found as 



(41) 



where the Lorentz factor 7 = l/vl — 1^ and the sign of tp is chosen to ensure that the velocity 76^ (y) of the 
electron fluid in the wave frame is non-positive. Thus, in the laboratory frame the speed of the electrons is 
slower than the phase speed of the wave, except at wave-breaking where the electrons catch the wave. 

In the present analysis we assume that the electromagnetic field is due entirely to the electron fluid and 
ion background. The magnetic field vanishes and the only non-zero component of the electric field is in the 
direction; it follows that the Maxwell 2-form F is 



F = Edx° A dx^ 

where E is the x^ component of the electric field. 

Since V and F are elements of the subspace of forms on Ai generated by {dx'^, dx^} it follows 



VA(dV- —F) 
m 







and, using 



with g{V, V) 



1, we obtain 



dV=^F 
m 



and (gniliailll) yield 
^ 1 m d/i 

Since dC A F ^ and di^ = it follows 
OCem 



d 



dY 



F = 



and inserting (gSl HO]) into @ yields 



dC V dX 



dfA 
Y=odC, J 
where 

r = 

follows from ^ and (HH) has been used. Hence 



- 7^ 



d 

dC 



,d&EM_ diJ 

' dX y^,dC 

The product of and dii/d( may be written 
d 



- 1 



dC 



5£em 
^^^^ 



EM 



mZnion{v\/ - 72 - n) 



Y=0 



usmg 

X = E^, 
which follows from 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 



(48) 



(49) 



(50) 



(51) 
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I 



III c 



II 

Figure 1: The dashed curve shows ^ versus C, and the soUd curve shows E versus Q (not to scale). Points of 
intersection of E with the ^-axis are labeled I and III, and II is a turning point of E. Using (|45|) and g < it 
follows dfi/dC, and E are of opposite sign. 



3.1 Maximum amplitude oscillation 

The square root in the right-hand side of (|49p places a lower bound on Suppose that over an oscillation the 
component ji oiV attains its lowest possible value /xi = 7 at ^ = Ci (see figure [T|) . Since /ii is a turning point of 
/i it follows E{C,i) = 0. At C = Cn the electric field i?(Cii) — ~Eraax where i?max is the amplitude of the largest 
possible oscillation (see figure [T|). Integration of (|50|) over [Ci,Cii] yields 



9£em ^ ^ 
oX 



Y=0 Y=G 



= mZ7iion(7 - 1) (52) 



For a standard cold Maxwell plasma £em — X/2, which inserted into ([5^ yields the well-known wave- 
breaking limit due to Akhiezer and Polovin |10| 



E^l^ = v/2mZnion(7-l) (53) 

mUJnC 



V2(7-l) (54) 



191 

for a multiply-ionized relativistic cold plasma, where 



= \ (55) 

is the plasma frequency and the speed of light c and permittivity of the vacuum eg have been restored. 
Using £em = -Cbi and ([TBI EH it follows 

i?,nL = - Jl-[^^iE^L)y2 + l]-^ (56) 

for a cold Born-Infcld plasma and limK-1.0 -^max = -^max as expected. 

The choice k — Eo^'o/I'J'I rnade by Born and Infeld, where rg is the classical radius of the electron, leads 
to K ^ 10~^^ni/V. However, quantum electrodynamic vacuum effects are expected to become important 
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at electric field strengths ~ lO^^V/m, and if Born-Infeld theory plays a role at the classical level it follows 
-^max ^ 10^^ V/m and KiJ^jj,^^ ^ 1. The next generation of laser systems (such as Extreme Light Infrastructure 
(ELI) [U and High Power Laser Energy Research system (HiPER) [IB]), are expected to offer intensities 
^ 10^5 W/cm^ corresponding to electric field strengths ~ lO^^V/m. 

Although in practice kE^^^ ^ 1 it is worth noting that although the limit of a.s v ^ c docs not exist, 

using it follows 

]iTnEZ. = - (57) 

•u— >-c K 

which is the value of the electric field of a static point electron evaluated at the location of the point electron 
in Born-Infcld electrodynamics [3]. 

3.2 Period and frequency of mciximum amplitude oscillations 

The period A of the maximum amplitude oscillation of a cold Born-Infeld plasma is obtained from the solution 
to (j49|) with £em = ^bi and the initial conditions ^(Ci) = 7 and c'yu/rfClc=Ci ~ '^h'^ particular first integral 
of (j49p satisfying the initial conditions on ^ may be written 

= 5^ (l - [^'rnZn-Uv^J^^ - + l) + 1]"') (58) 

and. using (|58p . consideration of the stationary points of fi yields 7 < /i < 7'^(1 + u^)- Furthermore, since 
/^(Ci) = 7 it follows ^ > for < ^ < Cm where n{Ciu) = 7'^(1 + v'^) is the maximum value of /i (see figure[T|). 
Thus, using (j58p the period A of the maximum amplitude oscillation is 



A = 2(Ciii - Ci) 

7-'(l+i'^) 



k — dfi 



(59) 



where k = Kmcjjp/\q\. 

A lower bound on A may be determined by noting < (uyT^^ — 7^— + 7) < 7 — 1 for 7 < < 7^(1 + f ^). 
It follows 

^ > ^ , (60) 



- (t;^];^^ - + 7) + l]-2 ['*'(7-l) + l] 

and thus 



7^(1+1''') 

2 



A > / k — du 
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Hence, using (|57l l6T|) it follows A diverges as least as fast as 7 in the limit w c. Thus, although the 
Lorentz force on an electron trapped in a maximum amplitude wave is finite, the estimate Arwyv^ /\q\ of the 
work done on the electron over a wave period diverges in the limit v ^ c and, in principle, it is possible to 
accelerate electrons to arbitrarily high energies in the present classical theory. Furthermore, consideration of the 
electromagnetic stress-energy-momentum tensor 7em arising from metric variations of the first term in psp [3] 
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yields an electromagnetic mass-energy density qem (in the lab frame) that diverges in the limit w — > c for a 
wave on the verge of breaking. The tensor 7em satisfies 



Tem{T, U) = ^^^T [iuF A*F-FAlu*F 

where T, U are arbitrary vectors, and the 0-form qem is 

£'EM = TEM{d/dx°,d/dx"). 
Using dm El [Ml) it follows 



Y=0 



dX 



dY 



giT, u) 



and the Born-Infeld energy density q-qi obtained using (|64l [T6| 
1 1 



k2 VI - 



(1- 



(62) 



(63) 



(64) 



(65) 



diverges in the limit £' -> 1/k. 

The behaviour of (1591) for ^ 1 and 7 3> 1 may be extracted by changing the integration variable to 
X — m/t^- Using ([5^)) it follows 



A = — 



27 



4 



l + -/i2/ + C'(^^/ 

4 



(66) 



where 



/ = 73(wv/x' -7-4 - X + 7"') 
= 7(1-^%) +0(7"')- 



(67) 



Since ^ = — wa;*^, the (angular) frequency wbi of the electrostatic oscillations measured in the lab frame is 

(68) 

and (ElETlEg) yield 



2'KV 27r 

^BI = — = — +0(7"') 



WbI ~ t^AP 



1 - 



2q 
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(69) 



where the speed of light c has been restored and wap is the (angular) frequency of electrostatic oscillations of 
a cold Maxwell plasma for 7 1 [10] , 



WAP 



2^27 



(70) 



4 Conclusion 

In practice, a typical laser-plasma experiment currently has 7 ^ 10 — 100 and cjp ~ 10^* rad/s; using k 10^^^ it 
follows (wBi/wAp — 1) 10"^"^ — 10^^2. Although this value is currently out of reach of terrestrial experiments, we 
suggest that it may be possible to investigate the ramifications of ((551 [CTl [S5)) in the environment of astrophysical 
objects such as quasars and magnetars. In any case, a number of significant items have been omitted from the 
above calculation, such as the finite temperature of the plasma and ion motion. Future work will include more 
thorough evaluation of the consequences of non-zero k in the context of more comprehensive plasma models. 
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